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We evaluated homologous recombination deficient (HRD) phenotypes in epithelial ovarian cancer (EOC) 
considering BRCAl, BRCA2, and RADSlCm a large well- annotated patient set. We evaluated EOC patients 
for germline deleterious mutations (n — 899), somatic mutations (n — 279) and epigenetic alterations (n = 
482) in these genes using NGS and genome-wide methylation arrays. Deleterious germline mutations were 
identified in 32 (3.6%) patients for BRCAl, in 28 (3.1%) for BRCA2 and in 26 (2.9%) for RAD51C. Ten 
somatically sequenced patients had deleterious alterations, six (2.1%) in BRCAl and four (1.4%) in BRCA2. 
Fifty two patients (10.8%) had methylated BRCAl or RAD51C. HRD patients with germline or somatic 
alterations in any gene were more likely to be high grade serous, have an earlier diagnosis age and have 
ovarian and/or breast cancer family history. The HRD phenotype was most common in high grade serous 
EOC. Identification of EOC patients with an HRD phenotype may help tailor specific therapies. 



Ovarian cancer represents the fifth most common cause of cancer mortality in women due in part to the 
advanced stage at which patients typically present, with an estimated 14,030 deaths and 22,240 cases in 
2013 in the United States'. The most common ovarian cancer is invasive epithelial ovarian cancer, with 
five common histological subtypes: high grade serous (70%), endometrioid (10%), clear cell (10%), low grade 
serous (5%), and mucinous (3%). A family history of breast or ovarian cancer in first degree relatives increases the 
risk for ovarian cancer"^; family-based studies have revealed high- and moderate-penetrance genes including 
BRCA1\ BRCA2\ DNA mismatch repair genes=, RAD51C '', RADSID" and BRIP1% and case-control studies have 
identified eleven common variants associating with modestly increased risks' ' \ Estimates of the contribution of 
germline BRCAl and BRCA2 mutations to EOC vary widely from 5% to 20%"'; somatic mutations in BRCAl and 
BRCA2 occur less frequently (in 2%-8% patients)" '" and, as has BRCAl methylation" "'^ have been reported in 
ovarian cancer patients with no family history^^ 

Ovarian cancer patients who are BRCAl and BRCA2 germline mutation carriers have been reported to have an 
improved outcome compared to non-carriers and the mechanism underlying this benefit has been hypothesized 
as a high response rate to platinum agents, particularly among patients with high grade serous histology^"*"^'. This 
is because patients with germline BRCAl or BRCA2 mutations have an impaired ability to repair double stranded 
DNA breaks through homologous recombination'^'*. More recently, favorable responses to poly (adenosine 
diphosphate-ribose) polymerase (PARP) inhibitors have been noted in ovarian cancer patients carrying dele- 
terious germline BRCAl and BRCA2 mutations^'. However, some ovarian cancer patients without germline 
BRCAl and BRCA2 mutations also respond to PARP inhibitors''", suggesting that broader dysfunction of homo- 
logous recombination than BRCAl and BRCA2, a homologous recombination deficient (HRD) phenotype may 
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be important. HRD may be the result of deleterious germline muta- 
tion in other homologous recombination genes such as RAD51B, 
RAD51C, or RAD51D or somatic mutation or silencing of the same 
homologous recombination genes'" ''^ 

We evaluated the prevalence of HRD in a large EOC case series 
and correlated this with clinical outcomes. HRD was determined by 
sequencing BRCAl, BRCA2 and RAD51C in germline and tumor 
DNA and assessing DNA methylation of cytosines at CpG sites 
surrounding these genes. 

Results 

Table 1 shows the characteristics of patients and Figure 1 shows the 
number of samples with each data type, after QC exclusions, 
included in this analysis. 



Germline Somatic 
Sequence Sequence 
(N=899) N=584 N=3 N=0 (N=899) 




N=13 



Frequency of methylation, germline and somatic mutations. Of 

the 899 germline-tested patients, deleterious mutations were seen in 
83 (9%): 32 (3.5%) in BRCAl, 28 in BRCA2 (3%), and 26 in RAD51C 
(2.9%; Table 2). Three cases carried mutations in more than one 
gene: two patients carried deleterious mutations in BRCAl and 
RAD51C, and one carried deleterious mutations in BRCA2 and 
RAD51C. For the 619 HGS cases, deleterious mutations were seen 
in 69 (11%), 27 (4%) in BRCAl, 24 (4%) in BRCA2 and 19 (3%) in 
RAD51C. One HGS patient carried a mutation in both BRCA2 and 



Tumor 
l\/1ethyiation 
(N=899) 

Figure 1 | Venn diagram representing data availability for N = 1063 
invasive EOC patients. 

RAD51C. Of the six Ashkenazi Jewish patients, one carried a 
germline BRCAl mutation (C61G) and none carried either somatic 
mutations or methylation. Table SI and S2 provides a complete 
mutation listing, after exclusions (see Methods section). Another 
187 women carried germline missense mutations, none of which 
had evidence for pathogenicity in the key clinical databases, 65 of 
these alterations were in BRCAl, 102 in BRCA2 and 32 in RAD51C. 
Twelve patients had a missense mutation in multiple genes. These 
variants were considered variants of unknown significance (VUS). In 
the patients with deleterious mutations, 13 patients also had a 
missense VUS that was identified. The remaining 629 patients had 
no deleterious or missense VUS mutations in BRCAl, BRCA2 or 
RAD51C. 

Two hundred and seventy nine patients with available tumor tis- 
sue were screened for somatic mutations. Ten tumors (4%) were 
found to have deleterious somatic mutations, six in BRCAl, and four 
in BRCA2 (Table 2). Table S2 provides a complete listing for the 
somatic mutations. Missense VUS were also detected with one in 
BRCAl and one in BRCA2. Two hundred and sixty seven patients 
had no deleterious or missense VUS mutations detected. No dele- 
terious somatic variations were identified in RAD51C. 

Forty five of 482 tumors (9%) exhibited a methylated BRCAl 
phenotype, and seven tumors (1%) showed RAD51C methylation 
(Table 2). BRCAl methylation was not detected in any individual 
carrying a germline BRCAl deleterious mutation; similarly, RAD51C 
methylation was exclusive of germline RAD51C mutation carrier 
status. However, there were two patients with a mutation and methy- 
lation of another gene; one patient had a somatic mutation in 
BRCA2 and BRCAl methylation and another carried a germline 
RAD51C mutation and had methylated BRCAl. 

Clinical and pathological analyses. In combination, 143 patients 
exhibited a HRD phenotype (82 due to germline mutation only, 
nine due to somatic mutation only, 50 due to tumor methylation 
only, one due to germline mutation and methylation, and one 
additional due to somatic mutation and methylation. Table 2). A 
total of 79 HRD patients showed alterations in BRCAl, 30 in 
BRCA2, 29 in RAD51C and five patients showed alterations in 
more than one gene. 

We examined the clinical features of HRD compared to patients 
with no abnormality detected (NAD) for those with both germline 
and somatic data (Tables 3, S3). HRD patients had younger age at 
diagnosis (P = 0.0001) and were more likely to be high grade serous 
(1 18/264, P = 0.0004) and have a family history of breast or ovarian 



Table 1 | Clinical characteristics of epithelial ovarian cancer 
patients (N = 1 063) 



N(%) 



Mean age at diagnosis, years (SD) 


62.4(11.7) 


Histology 




High-grade serous 


735 (70%) 


High grade endometrioid 


73 (7%) 


Low grade endometrioid 


67 (6%) 


Clear Cell 


69 (7%) 


Low-grade serous 


34 (3%) 


Mucinous 


29 (3%) 


Other/Unknown 


56 


FIGO Stage 






158 (15%) 


II 


62 (6%) 


III 


637(61%) 


IV 


194 (19%) 


Unknown 


12 


Debulking category 




Optimal 


848 (85%) 


Sub-optimal 


150(15%) 


Unknown 


65 


Grade 




Low grade 


143 (14%) 


High grade 


896 (86%) 


Unknown 


24 


Menopause status 




Pre-menopausal/peri-menopausal 


1 69 (20%) 


Post-menopausal 


672 (80%) 


Unknown 


222 


Ascites present at surgery 




Yes 


524 (63%) 


No 


307 (37%) 


Unknown 


232 


Ashkenazi ancestry 




Yes 


6 (1%) 


No 


498 (99%) 


Unknown 


559 



Median survival, years (range) 2.9 (0.01-10) 

Median follow-up time among living, years (range) 4.5 (0.01-1 0) 

Deaths 638 (60%) 

N (%), unless otherwise specified. 
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Table 2 | Summary of germline mutations, somatic mutations. 


and methylation 


of BRCAl, 


BRCA2 an 


d RAD51C 










Gene 








N Tested 


BRCAl 


BRCA2 




Combined HRD 


Oermline delererious mutation 


poo 
oyy 




z o 


zo 


Oo 


Somatic deleterious mutation 


ovo 


o 


4 


rv 
U 


1 u 


Methylated 


482 


45 




7 


52 


No deleterious mutation, not methylated 


NA 


237 


261 


263 


NA 


No deleterious germline mutation, unknown somatic mutation 


NA 


592 


606 


604 


NA 


and/or tumor methylation status 












Not methylated, germline and somatic mutation status unknown 


NA 


151 


164 


163 


NA 


Total 




1063 


1063 


1063 


143 



Deleterious mutations were frame shift insertion/deletion, splice site, rare missense supported by multiple metfiodsto be damaging; NA, not applicable; one case carried germline deleterious mutations for 
both 6RO\2and RAD51C, and two patients carried germline deleterious mutations for BRCA} and RAD51C, thus the combined number of germline deleterious mutation carriers is less than the sum of 
mutation carriers for each gene; one case carried a RAD51 C germline deleterious mutation and methylated for BRCAl and another case carried a BRCA2 somatic deleterious mutation and methylated for 
BRCAl , thus the total number of patients with a HRD phenofype is less than the sum of germline mutation carriers and patients with a somatically mutated or methylated gene. 



cancer (P = 0.002) than patients with NAD. Interestingly, patients 
with high grade (10/27) endometrioid EOC were HRD, compared to 
2/23 patients with low grade endometrioid EOC. One of the most 
common pathology challenges is whether high grade endometrioid 
are in fact high grade serous EOC''\ Other covariates (debulking, 
peritoneal cytology, the presence of ascites and smoking and alcohol 
use history) were not significantly different between the two groups. 
We also examined clinical features within the HRD patients and 
found only family history of breast or ovarian cancer was more 
common in mutated than methylated patients (P = 0.0001, data 
not shown). 

Survival analyses were carried out on those patients with HRD and 
those with NAD for which all data (germline, somatic mutation and 
methylation) were available (n = 356). Forty percent (143/356) of all 
the patients were HRD, while 44.7% (1 18/264) of HGS patients were 
HRD. Overall survival for all patients was associated with the HRD 
phenotype (adjusted HR 0.73, 95% CI 0.54-0.98, P = 0.04, Table 4), 
but not within each of the defective classes, genes that were mutated 
at the germline or somatic level or were somatically methylated 
(respectively, adjusted HR 0.81, 95% CI 0.57-1.15; adjusted HR 
0.56, 95% CI 0.20-1.54, adjusted HR 0.65, 95% CI 0.42-1.02, P = 
0.18, Table 4), or when outcomes were examined by germline muta- 
tion carrier status (P = 0.09, Table 4). The hazard ratios data are 



similar to those previously reported in a large germline study^*". 
Results were similar whether VUS were included in the NAD cat- 
egory or excluded (data not shown). Analysis of time to recurrence 
showed no association with any of the HRD phenotypes (Table S3). 
Survival analyses were also conducted restricted to high grade serous 
(n = 618 for germline carrier analysis, 264 for HRD analysis) and 
showed an overall survival differences between HRD and NAD 
(adjusted HR 0.70, 95% CI 0.51-0.97, P = 0.03, Table 4, Figure 2), 
but not between germline mutation carriers and-non carriers 
(adjusted HR 1.07, 95% CI 0.65-1.76; adjusted HR 0.76, 95% CI 
0.45-1.30, adjusted HR 0.66, 95% CI 0.37-1.18, P = 0.40, Table 4, 
Figure 2). Survival at 5 years was similar (adjusted HR 0.68 95% CI 
0.48, 0.97, P = 0.04, data not shown). Analysis of time to recurrence 
was null (Table S4, Figure S3). Finally, we examined HRD pheno- 
types by debulking status in germline mutation carriers, as Alsop 
et reported that improved survival was not evident in patients 
that were sub-optimally debulked. However, we observed no differ- 
ence in germline mutation carriers stratified by debulking status 
(data not shown). 

Discussion 

Using a large homogeneous collection of ovarian cancer patients 
with clinical follow up unselected for family history, we evaluated 



Table 3 | Comparison of the clinical characteristics of HRD to NAD patients 

HRD(n = 143) NAD (n = 21 3) Rvalue 



Age at Diagnosis, years 



Mean (SD) 


58.4 (10.6) 


63.1 (1 1.5) 


Histology 






High grade serous 


1 1 8 (45%) 


146 (55%) 


High grade endometrioid 


1 0 (37%) 


1 7 (63%) 


Low grade endometrioid 


2 (9%) 


21 (91%) 


Clear Cell 


5 (28%) 


1 3 (72%) 


Low grade Serous 


2 (15%) 


1 1 (85%) 


Mucinous 


0 (0%) 


5 (100%) 


Other/Unknown 


6 (100%) 


0 (0%) 


Grade 






High 


1 35 (44%) 


1 74 (66%) 


Low 


5 (11%) 


39 (89%) 


Unknown 


3 (100%) 


0 (0%) 


First Degree Family History of Breast or Ovarian Cancer 






No 


104 (36%) 


1 83 (64%) 


Yes 


39 (57%) 


30 (43%) 


First Degree Family History of Ovarian Cancer 






No 


132 (39%) 


207 (61%) 


Yes 


1 1 (65%) 


6 (35%) 



0.0001 
0.0004 



<0.0001 

0.0020 
0.0344 



Only patients with known HRD status (germline mutation, somatic mutation and somatic methylation) available were included; VUS carriers are included in NAD patient group; stage, debulking, presence of 
ascites, peritoneal cytology, smoking history, enrollment year, alcohol use were not significantly different between HRD and NAD patients. 



SCIENTIFIC REPORTS | 4 : 4026 | DOI: 1 0. 1 038/srep04026 



3 



-a 

D 

o 



o 



g 
b 



Q. 
O 



Q. 
Q 



_Q 



D 
c 

< 



o 



o 



o 
0 



X 



o 
0 

g 

z 





K 
a 






CN 

q 




0 0 00 

K CO 




0^ 














>J~) 




•d CO 


CO 




ut" K 
0 CO 




d 


0 


do. 


d 




d 0,0 




0 

K 




0 CN 

00 


CO 

<) 




K 0 0 

q K <i 




d 




d d 


d 




— d d 


o 




K 

CN 






0 
CN 




d 




d 






d 






0 
q 




ijT"co" 






cn"?To" 

CN CO — ; 














0' ut" 
^ ^ CO 




d 


"0 


do. 


d 


"0 


2.2.2. 


K 


^ K 

00 




K CN 
K K 0 




d 




d d 


d 




d d d 




0^ 


0^ 






^ 
0^ 


^ 

0^ 0^ 0^ 


CO 

0 


<5 


00 


<J~) CO 
0 ^ 




K 
•0 


CO — l\ 

000 


0 ^ 

0 K 


0 ^ CO 


K 
CN 


K K ^ CN 
0 ^ ^ ^ 
CO 




00 


0 00 K 

^ 0 




0 


K CO 00 
CM CN ^ 


0 




00 






a 
0 




CD 




d 






d 






00" 
a 




IT) ^ 


cn" 
q 




uTo"^ 

00 — 0 




d 


















t\ o" 

CM 


CN 




CN 00" C0~ 
K CO CO 


"0 


d 


"0 


do. 




"0 


0,0,0 




K 




^ <5 
00 lO 






0 00 

^ 0 




d 




d d 


d 




^ d d 


K 
CM 




K 






K 

CO 




CD 




d 






d 










^cn" 

CN CO 


0" 

CO 




^^cn" 

lO CO 




•0 




ut" 

0 CN 


d" 




CN o" 

0 CO 


"0 


d 


"0 


do. 


d 


^0 


d 0,0 


00 


0 K 
00 <3 


00 
00 


0 CS CO 
C> K 0 




d 




d d 


d 




d d d 




0^ 










^ 


CN 




CN 


c> 


K CO OM CN 


CO 


00 


CO 


00 


CN 
CO 


0 C> ^ CN 
0 ^ ^ ^ 


n 

CN 


CO 


CO 
CN 


CN 0 

00 


0 


•0 0 K CO 
^ CO CM CN 
00 



5 E 



.0 
D 

B J 

0 



0 s 



2 E 



* .E 

Q) 



t£ a 



Z X I 



D ~ CN 

V < <r JO 
a; 0 gjf?; o^yS 



cn Z 



5 a 



the clinical and pathological characteristics associated with a variety 
of HRD phenotypes. Deleterious germline mutations in BRCAl, 
BRCA2 or RAD51C were seen 9% of the subjects, BRCAl or BRCA2 
in 6.7% which is within the ranges of those previously reported (5.8%- 
13.2%)"". Somatic mutations were less common (3.6%), but 1 1% of the 
tumors analyzed exhibited BRCAl or RAD51C methylation, compar- 
able to the prevalence observed in 319 high grade serous ovarian 
cancers analyzed by the Cancer Genome Adas". In total, 143 of the 
patients had an HRD phenotype. A common phenotype has been 
reported for sporadic cancers that share the characteristics of tumors 
associated with germline BRCAl and BRCA2 mutations''''^ and 
includes RAD51C in which germline mutations have been identified 
in patients of familial breast and ovarian cancer''"'"'', as has methyla- 
tion'""". Additionally, Abkevich et al*" report an LOH-based score for 
HRD which included promoter methylation of RAD51C. Here, we 
analyzed patients regardless of family history, included germline as 
well as tumor testing for more comprehensive analysis and defined 
patients as HRD if any deleterious alterations in BRCAl, BRCA2 and 
RAD51C genes were found. 

An HRD phenotype was seen in 45% of high grade serous EOC 
(69% of patients had HQS EOC), with a smaller proportion 
observed in low grade serous and the other histologies. This is 
consistent with observations of others'''^*' '''. While we did observe 
a survival advantage for those patients with an HRD phenotype, 
analysis did not reveal any overall survival differences between 
BRCAl or BRCA2 mutation carriers and non-carriers, unlike prior 
studies. This was unexpected given the literature suggesting longer 
progression free survival or overall survival for EOC in mutation 
carriers; however, this advantage appears to be short term. Most 
studies were out to 60 months only and McLaughlin et al*" showed 
that an early (up to five years) survival benefit was not evident 
beyond eight years'". Similarly we do not show a long term sur- 
vival advantage, thus follow up time may explain the apparent 
discordance. It may also be that the number of germline BRCAl 
and BRCA2 mutations carriers was too small and thus underpow- 
ered. Additionally, while care was taken to classify the mutations 
according to the best available evidence, the power of our analysis 
will have been reduced by any misclassification of HRD. While the 
TCGA concluded that methylated patients more closely match 
survival in wild type BRCA patients, we observed that methylated 
cases may have a survival benefit similar to that of the germline 
mutation patients. All categories, while not reaching statistical 
significance, showed a trend towards a survival benefit; had 
methylation not been thus, then combining with mutation positive 
patients would have biased the testing procedure towards the null. 
Biologically, BRCAl inactivated by methylation would result in a 
similar defect as a deleteriously mutated allele. Our study involves 
patients at one site, and all data were similarly processed while the 
TCGA study involves samples from multiple centers with multiple 
data sets, which may also account for the differences. 

The high frequency of HRD defects in high grade serous EOC 
suggest screening for this phenotype should be considered in clinical 
care, particularly for new therapies, such as the PARP inhibitors that 
may be most efficacious within this phenotype. 

Strengths of this report include detailed methylation phenotyping 
of a subset of patients and sequencing of tumor DNA. Limitations 
include lack of methylation and somatic mutation data on up to 67% 
of patients, not confirming the detected mutations in every patient 
and inclusion of only three of more than ten genes now considered to 
contribute to the HRD phenotype''''. It is also critical to note that most 
of these patients were treated with the current standard of care at 
time of diagnosis (carboplatin and/or taxane); upcoming studies will 
focus on those treated with PARP inhibitors. Other future work 
includes sequencing of additional DNA repair genes in order to more 
comprehensively identify those patients with defective homologous 
recombination pathway abnormalities. 
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Figure 2 | Overall survival for high grade serous EOC by HRD and type of alteration. (A). Red: HRD; Black: NAD. (B). Above is a depiction of samples 
with type of alteration in each gene; below are the Kaplan Meier curves. Red: germline deleterious mutation; Green: somatic deleterious mutation; Blue: 
somatic methylation; Black: no germline mutation, somatic mutation or methylation. 



Methods 

Patient population and biospecimens. Women witli EOC {n — 1063) provided 
written informed consent between 1992 and 2011 using IRB approved processes at 
the Mayo Clinic in Rochester, MN. Use of patient derived samples in this study was 
also approved and the work described was carried out in accordance with these 
approved guidelines. Eligible participants were women aged 20 years or older with 
diagnosis of pathologically confirmed invasive primary epithelial ovarian, fallopian 
tube or primary peritoneal cancer. Patients were followed for vital status using 
electronic medical records and linkage to vital statistics sources; the median time of 
follow up was 4.5 years, with a range of 0. 1- 1 0 years. Most patients were treated with a 
platinum agent (carboplatin or cisplatin) and taxane, the standard of care at time of 
diagnosis. Diagnosis (histology and grade) were reviewed by a gynecologic 
pathologist. Peripheral blood (n — 899) was used a source of germline DNA'*''. Pre- 
chemotherapy tumor tissue (up to n — 482) was snap frozen immediately following 
surgery and stored at — 80''C; tumor DNA was prepared from cryostat sections using 
Qiagen PureGene chemistry (Qiagen Inc, Valencia, CA, USA) and quantitated using 
Trinean (Trinean, Gentbrugge, Belgium). Tumor RNA was isolated from fresh frozen 
samples using the Qiagen RNEasy protocol (Qiagen Inc, Valencia, CA, USA) and 
quantitated using a Nanodrop Spectrophotomer (Agilent Technologies, Santa Clara, 
CA, USA). Somatic methylation analyses preceded the mutation screening, thus 
samples with insufficient DNA for both analyses were not sequenced. 

DNA sequencing and bioinformatics. Germhne DNA (n = 899, 1.25 \i\ DNA at a 
concentration 75 ng/jil) and tumor DNA (n — 279) were sequenced at the University 
of Cambridge, Cambridge, UK at BRCAl, BRCA2 and RAD51C. The procedure 
entailed design of PCR primers, testing and optimization using Eluidigm Access 
Array (Eluidigm Corp., San Erancisco, CA, USA). Briefly, the method used a 4-primer 
amplicon-tagging scheme. Tagged, target-specific primer pairs were combined with 
sample-specific primer pairs that contained barcoding sequences and the adaptor 



sequences used by the lUumina sequencing system. Using sample-specific barcodes, 
all PCR products generated in the 48.48 Access Array Integrated Eluidic Circuit were 
unique and pooled together to run in a single sequencing experiment. Sequencing 
reactions were carried out on an lUumina GAIIX next-generation sequencer 
(lUumina, Inc., San Diego, CA, USA). 

GENOME_GPS 1.0 was used for analysis of DNA sequence data including read 
alignment, variant detection, and annotation (Figure SI). Reads were aligned to 
genome build hgl9 using Novoalign and realignment, recalibration used the Genome 
Analysis Toolkit (GATK) v. 1 .6.7^. The percent of reads mapped on target was used to 
exclude seven germline samples with <80% of amplicons covered at 20X and twenty 
one tumor samples with <80% of the amplicons covered at 20X. To avoid false- 
negative deleterious carrier status assignments for the somatic mutation data, non- 
carrier patients were excluded if either germline or tumor sequencing data exhibited 
<80% on-target sequence coverage at 20 X or greater [mean coverage, and B""'^ 
quartiles: RAD51C 138.76 (121.68, 150.94); BRCAl 164.81 (142.59, 179.47); BRCA2 
157.65 (133.68, 174.13); all three genes, 158.20 (136.06, 173.04)]. 

Germline variants were called through GATK Unified Genotyper. Germline filtering 
used genotype quality greater than 20 and read depth greater than 20. Somatic single 
nucleotide variants were called using SomaticSniper'*'\ whereas insertions and deletions 
were called by GATK Somatic Indel Detector''"*. Somatic variants were considered with 
a depth of at least 20 X coverage for the germline and tumor DNA at the given site as 
well as a quahty score ^ 20 and excluded if close to indel calls, aligned to multiple 
positions and/or occurred in a known repeat region. Indel calls required ^ 10 sup- 
porting tumor sequencing reads, and no reads supporting the indel in the germline. 

Each variant was annotated using the Targeted RE-sequencing Annotation Tool 
(TREAT)'"'. Missense coding variants were functionally annotated by snpEFE (http:// 
snpeffsourceforge.net/) and PolyPhen-2 to predict biological effects. Somatic 
variants were considered deleterious if they were stop-gains (nonsense) or occurred at 
splice sites junctions, along with frameshifts and insertion/deletions as deleterious. In 
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addition, missense mutations were hand-curated using LOVD (www.LOVD.nl), 
lARC (www.iarc.fr), HGMD (www.hgmd.cf.ac.uk), and Alamut (littp://www. 
interactive-biosoftware.com/ software/alamut/features) databases. Variants for which 
deleteriousness could not be determined with confidence, or which were known 
polymorphisms, were considered VUS. 

DNA methylation. Tumor methylation status for 482 patients was examined using the 
lUmnina Infinium HumanMethylation450 BeadChip as previously described^^. Positive 
bisulfite modified (BSM) controls {Sssl treated DNA, CpGenome Universal Methylated 
DNA, Millipore, Billerica, MA, USA, and placental DNA), negative BSM controls 
(unmodified Sssl treated DNA and WGA DNA), and a CEPH DNA (Coriell Institute), 
were included in each 96 well plate. Validation of nine CpG loci, including one at 
BRCAl (cg25288140), was accomplished by Pyrosequencing, as previously described^' 
and su^ested that BeadChip data were satisfactory; Pearson correlation of methylation 
values between BeadChip and Pyrosequencing assays ranged from 0.57 to 0.96. 

To dichotomize tumor methylation status for each gene, we used expression array 
data from 171 patients to identify CpG probes which negatively correlated with 
expression. Expression data were generated using 750 ng of high quality total RNA 
and Agilent Whole Human Genome 4 X 44K Expression Arrays (using a mixed 
reference of 107 tumors) as previously described^^. We estimated pairwise Spearman 
correlations between gene expression and CpG methylation probes within 20 kb of 
each gene to identify CpG probes most associated with expression. For each probe, we 
confirmed that each CpG probe was close to known promoter and active regulatory 
elements defined by ENCODE project using the UCSC Genome Browser (Figure 
S2A). Thresholds for dichotomization of methylation status were defmed based upon 
the median methylation value of the selected probe set. For example. Figure S2B plots 
BRCAl expression versus the median of N — 21 BRCAl CpG probes; based on this, 
patients with mean BRCAl CpG methylation > 0.15 were considered BRCAl 
methylated. A similar process was used for six CpG probes in RAD51C; however, no 
BRCA2 CpG was clearly associated with altered gene expression, thus BRCA2 was not 
part of methylation analyses. 

Association analysis. HRD patients were defmed as those with deleterious 
germline mutations, deleterious somatic mutations or methylation of BRCAl, 
BRCA2, or RAD51C. Clinical characteristics of HRD and patients with NAD were 
compared using Fisher exact tests. Survival analysis with HRD status was 
conducted using multivariable Cox proportional hazards modeling and, for 
adjusted analyses, included age, stage, debulking, menopausal status, grade, and 
ascites, estimating hazard ratios (HRs) and 95% confidence intervals (CIs). Time 
at risk was defined from the date of diagnosis to date of death allowing for left 
truncation with right censoring at last date of follow up to 10 years. In addition to 
overall survival, we considered time-to-recurrence or death from any cause. We 
analyzed survival using the HRD phenotype defined by germline status only and 
also by combined germline, somatic, and methylation status (using the subset 
with complete data). Analyses were for all patients and among patients with high 
grade serous EOC and were conducted both with VUS considered to be non- 
carriers or excluded. Each analysis used all available data (thus sample sizes vary) 
and was implemented in SAS, with corresponding Kaplan-Meier curves 
generated in R. 
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